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We examine vortex dynamics and pinning in layered superconductors using three-dimensional molecular dy- 
namics simulations of magnetically interacting pancake vortices. Our model treats the magnetic interactions of 
the pancakes exactly, with long-range logarithmic interactions both within and between planes. At the matching 
field the vortices are aligned with the pinning array. As a function of tilt angle for the pinning arrays a series of 
commensuration effects occur, seen as peaks in the critical current, due to pancakes finding a favorable alignment. 



In superconductors with periodic pinning ar- 
rays interesting commensurability effects occur 
when the periodicity of the vortex lattice matches 
the periodicity of the pinning lattice. Experi- 
ments |l],|| and simulations || so far have been 
done with thin film superconductors where the 
vortex lattice and pinning can be considered two- 
dimensional. The case of vortex lattices inter- 
acting with a periodic pinning array in a layered 
3D superconductor has not been studied. Such 
a system would correspond to an anisotropic su- 
perconductor such as BSCCO with a periodic ar- 
rangement of columnar defects. In this system 
the z-direction becomes important as the applied 
field or the pinning array is tilted. The dynamical 
effects of vortices moving in periodic pinning ar- 
rays in such a system have not been examined, in 
particular how the vortex lattice structure of the 
moving state differs from that of the pinned state. 
To study vortex pinning and dynamics in layered 
superconductors, we have developed a simulation 
containing the correct magnetic interactions be- 
tween pancakes [Q . This interaction is long range 
both in and between planes, and is treated using 
a rapidly converging summation method ||] . 
The overdamped equation of motion, for T = 0, 



i? p -I- id — Vj , where N v is the number of vortices 
and p and z are the distance between pancakes 
in cylindrical coordinates. The magnetic energy 
between pancakes is 



U(fH,j,z) 



where 
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for vortex i is given by ii = Ylj=i VU(pj 5 
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R = E x {x) = f~exp(-x/\)/f/dp? 

and eo = (47r^) 2 . The pinning is placed in 
a square array of parabolic traps with a radius 
r p much smaller than the distance between pins. 
The location of the pinning sites is the same in 
every layer corresponding to correlated defects. A 
driving force fd is slowly increased and the vor- 
tex velocities are measured. Here we consider the 
first matching field case where the number of vor- 
tices N v equals the number of pinning sites N p . 
We conduct a series of simulations in which the 
pinning sites are tilted at an increasing angle with 
respect to the z-axis. We will only consider driv- 
ing that produces vortex motion transverse to the 
direction of the tilt angle. We examine systems 
with 8 layers containing 64 vortices and pins in 
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each layer. Work for larger systems, varied fields 
and coupling strength will be presented elsewhere 

In Fig. 1(a) we present the critical depinning 
force /J as a function of tilt angle 8. Here /J p 
peaks at 8 = 0° when the pancakes are aligned 
with pins on all layers. As 8 is increased f% 
drops. For small tilt angles 8 < 5° the vortex 
lines tilt with the pins. For larger angles the vor- 
tex lines realign in the z direction. The depinning 
force fH will then remain low as only one pan- 
cake in the straight vortex line will be sitting at 
a pinning site. At 8 — 45° f% shows a peak of 
the same magnitude as the peak at 8 = 0. At 
this tilt angle, and also for any angle satisfying 
8 = tan -1 (ro) where n is an integer, the pinning 
sites are again aligned in the z-direction so that 
a vortex line can be formed that is also aligned in 
the z-direction with all the pancakes in a single 
vortex being able to sit in a pinning site. There 
are also peaks in f% at 8 = 26.6° and 56.3°. At 
these angles the pancakes again sit on all the pin- 
ning sites. The individual vortex lines now con- 
sists of half the number of pancakes as at 8 = 0.0°; 
however, there are now twice as many vortex lines 
with the pancakes from an individual vortex line 
being coupled in every other layer. The view from 
the z-direction as shown in Fig. 1 for these angles 
indicates that the vortex lattice is now rectan- 
gular with twice as many vortex lines as at the 
other angles. At 6 = 36.9° a smaller peak is ob- 
served. The vortex structure at this angle will be 
presented elsewhere Q. 

In (b) and (c) we show the vortex structures for 
the pinned phase and moving phase for 8 = 1.5° 
as seen from the z-direction. In (b) the vortices 
can be seen to stay aligned with the pins. In 
(c) for jd > fj p the vortices realign with the 
z-direction. Such a transition from a tilted to 
straight vortex lattice as a function of drive may 
be visible with neutron scattering experiments. 
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Figure 1. The critical depinning force /J p ver- 
sus the tilt angle 8 of the pinning sites. The vor- 
tex arrangements as seen from the z-direction are 
outlined for different tilt angles 8 = left and 
8 = 26.6 right, (b) shows the pinned vortex ar- 
rangement for 8 = 1.5° where the vortices stay 
aligned with the pins, (c) shows the moving vor- 
tex state for 8 = 1.5° where the vortices have 
realigned with the z direction. 
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